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Abstract: We show that the late time rolling of the Cubic Superstring Field Theory 
(CSSFT) non-local tachyon in the FRW Universe leads to a cosmic acceleration with a 
periodic crossing of the w = —1 barrier. An asymptotic solution for the tachyon and 
Hubble parameter by linearizing the non-local equations of motion is constructed explicitly. 
For a small Hubble parameter the period of oscillations is a number entirely defined by the 
parameters of the CSSFT action. 
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1. Introduction 

The combined analysis of the type la supernovae, galaxy clusters measurements and WMAP 
(Wilkinson Microwave Anisotropy Probe) data brings out clearly an evidence of the accel- 
erated expansion of the Universe @-|5|. The cosmological acceleration strongly indicates 
that the present day Universe is dominated by a smoothly distributed slowly varying Dark 
Energy (DE) . Recent results of WMAP || together with la supernovae data give a strong 
support that the present time DE state parameter is close to —1: 

w = -0.97^ 9 7 

or without an a priori assumption that the Universe is flat and together with large-scale 
structure and supernovae data w = — 1.06j^[J|. 

From a theoretical point of view the above mentioned domain of w covers three es- 
sentially different cases. The first case, w > —1, is achieved in quintessence models M, ||[ 
containing an extra light scalar field which is not in the Standard Model set of fields 0]. 



The second case, w = — 1, is the cosmological constant [|10|, |11|. The third case, w < — 1, 
is called a "phantom" one and can be realized by a scalar field with a ghost (phantom) 
kinetic term. In this case all natural energy conditions are violated and there are problems 
of an instability at classical and quantum levels [12, 13 1. 



Since experimental data do not contradict with a possibility w < — 1 and moreover a 



direct search strategy to test inequality w < — 1 has been proposed [14] a study of such 
models attracts a lot of attention. Some projects p5| explore whether w varies with the 
time or an exact constant. Varying w obviously corresponds to a dynamical model of DE 
(see [16] for a review) which generally speaking includes a scalar field 1 . 



1 Modified models of GR also generate an effective scalar field (see for example |l7| and refs. therein). 
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A possible way to evade the instability problem for models with w < —1 is to yield 
a phantom model as an effective one, arising from a more fundamental theory without a 
negative kinetic term. In this paper we develop in more details a cosmological SFT tachyon 
model |T^] . The model is based on an SFT formulation of a fermionic NSR string with the 
GSO— sector [19]. In this model a scalar field is the open string tachyon, which describes 



according to the Sen's conjecture p0|] a dynamical transition of a non-BPS D-brane 2 to a 



stable vacuum (see [22 1 for review). Since the concerned model is a string theory limit, all 
stability issues are related to a stability of a VSFT (a Vacuum String Field Theory, i.e. 
the SFT in a true vacuum) and one has to discuss only an application of this limit to a 



full string theory. There are general arguments [23] that there does not exist a local scalar 
field model for a phantom Universe without an UV pathology. In a recent paper 24] it has 
been proposed a phantom model without UV pathology in which a vector field is used. 

The scalar model we investigate in this paper is a nonlocal one. Our goal is a con- 
struction of an analytic solution to linearized Friedmann equations in this model at large 
times. A characteristic feature of this model in the flat background is a presence of a rolling 



tachyon solution [25, 26, 27]. This property persuades us to consider a fermionic string 



instead of bosonic one where such a solution does not exist [28, 29 1. However one might 
expect an existence of a rolling solution in bosonic strings in a non-flat case. The dynamics 
of a non-local tachyon scalar field on a cosmological background in the Hamilton-Jacobi 
formalism is studied in [30]. 



We find that during the late time evolution in the FRW Universe the tachyon goes 
to its minimum oscillating with an exponentially decreasing amplitude. This is similar to 
the flat case. Consequently the Hubble parameter goes to a constant, and the state and 
deceleration parameters go to —1 all oscillating around their asymptotic values. The DE 
state parameter w crosses the phantom divide w = — 1 during an evolution. 

Models with a crossing of the w = —1 barrier are also a subject of recent studies. 
Simplest models include two fields (one phantom and one usual field, see [£| and 

refs. therein). General ft-essence models |M| can have both w < — 1 and w ^ — 1 but 
a dynamical transition from the region w ^ — 1 to the region w < — 1 or vice versa is 



forbidden under general assumptions |35|] and is possible only under special conditions |36| . 

In our case a non-locality provides a crossing of the w = — 1 barrier in spite of the 
presence of only one scalar field. Hence the Universe exhibits an acceleration but because 
of oscillations quintessence and phantom phases change one each other with time. 

The paper is organized as follows. In Section 2 we setup the cosmological model 
which is an approximation of the CSSFT describing a non-BPS brane within the level 
truncation scheme in the FRW Universe. In Section 3 we present details of the tachyon 
dynamics in the flat case at large times where an approximation linear in fluctuations 
around a non-perturbative vacuum is valid. We compare this linear approximation with 
a numeric solution to full equations. In Section 4 we study the tachyon dynamics in the 
FRW background again using a linear approximation to the Friedmann equations and 
find out that the obtained solution describes an accelerating Universe. In Section 5 we 



2 DE models based on brane- world scenarios are presented in 
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discuss cosmological consequences of the obtained results and point out further directions 
of studying this type of models. 



2. Setup 

An action for the tachyon in the CSSFT [37, 38] in the flat background 3 when fields up to 
zero mass are taken into account is found to be p5[ 

1 f ( n' 1 p 2A 

Ssft = / dx u \ x ) - ^ V dMx)d v <P{x) + -<p\x) + — cj> 2 (x)u(x) 



g l o! J \ 2 4 3 

where 4>(x) is the tachyon field, u{x) is an auxiliary field, 

and A = — log ~ 0.2616. r/ is the flat Minkowskian metric, □ = rj^d^d^. For simplicity 
we will use hereafter a' = 1 units in which all fields, coordinates and the coupling constant 
g are dimensionless. 

An auxiliary field u(x) can be integrated out to yield 

Stach = ^ j dx (^- l -^d^(x)d u m + \<f{x) - Q {X)\ . 

A reasonable but an ad hoc assumption 4 that u has no the tilde simplifies the last term 



in this action. Namely, under this assumption and a rescaling x — ► 2y2Xx 1 (j) — > 
and <7 — > 12Ae~ 2A g D the action for the tachyon becomes 

Stach, approx = ^2 j dx (~f lf v d„(f>(x)d v <j)(x) + ^(x) - I( e I G 0) 4 (x)) (2.1) 

where £ 2 = Jr ~ 0.9556. The last term in this action contains an infinite number of 
derivatives. Just due to this nonlocal factor a novel behavior in a dynamics of the tachyon 
field appears [25]. 

Cosmological scenarios with our Universe to be considered as a D3-brane embedded in 
10-dimensional space-time was proposed in and a dynamics of this brane is given by 
the following covariant version of action in a non-flat space 

S = j efcv^^ + ^-^a^^arJ + i^CxJ-^xJ-A^ (2.2) 



where 



3 We always use the signature ( — ,+ ,+ ,+ ,...). 

4 As it was shown in |23, Eq] this assumption can be applied to study a special class of rolling solutions 



and related questions. 
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Here g is the metric, k is a gravitational coupling constant and we choose such units that 
it is dimensionless, A is a constant. 

In the present analysis we focus on the four dimensional Universe with the spatially 
flat FRW metric which can be written as 

g^y = diag(-l,a 2 ,a 2 ,a 2 ) 

with a = a(t) being a space homogeneous scale factor. In this particular case □ is expressed 
through a as 

1 

7 2 

where H = a/a is the Hubble parameter and the dot denotes the time derivative. 
The Friedmann equations for the space homogeneous tachyon field have the form 



D g = -d 2 - 3Hd t + -dl 



3# 2 = 4 ( ^ 2 -U 2 + l -^ + £ 1+ £ 2 + k), (2.3a) 



9 2 V 2 r 2 r 4 
9 2 o\ 2 



H=^(-^ 2 -£ 2 ) (2.3b) 



where 



Ex = -- f ds(el sV ^ 3 )Ve-l sV ^>, £ 2 = -- f ds{d t e\ sZ y <S> 3 )d t e-^ sT '$ 
8 Jo 8 Jo 

with 

$ = et v <f>, V = -<9f - 3H(t)d t . 
The equation of motion for the tachyon is 

(fV + l)e-^ v ^ = $ 3 . (2.4) 

The latter equation is in fact the continuity equation for the cosmic fluid. To see this 
explicitly the following operator relation is useful 

A J ds ([e sXU °n g <t))Pe {1 - s)XU ^ - (e sAD ^p 9 Pe (1 - s)An ^ 
=(e AD "0)PV - </>Pe XD ° ifi - A / ds{e sXU ^)Tie {l - s)XU ^ 



where P and n are operators satisfying [\D g , P] = ir. 

Equations ( |2.3| ) are complicated because of the presence of an infinite number of deriva- 
tives and a non-flat metric. Before we are going to present a construction of an asymptotic 
solution to these equations to gain an insight into the problem of dealing with an infinite 



number of derivatives we shortly review the results of an analysis of action (2.1) in the flat 



space 25, 26, 27j 
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3. A Late time Rolling Tachyon in the flat background 



Approximate action for the tachyon (^J) in the flat background was already studied in 
p5| , |26| , p7|] 5 and we will present the results relevant for the sequel. The equation of motion 
for space homogeneous configurations of the tachyon field is found to be 



(-fd? + 1) e^$>{t) = $(tf. (3.1) 



where <3? = e~» 9 ^ (ft. The operator e s<9 * for positive s can be rewritten in an integral form 
and this allows further to implement an iteration procedure which is capable to find a 
numeric solution to the problem. 

If we are interesting in the Rolling Tachyon solution then we have in mind the following 
picture. The tachyon field starts from the origin (may be with a non-zero velocity), rolls 
down to the minimum of the tachyon potential 6 and eventually stops in the minimum. In 
our notation the minima are located at <3?o = ±1. For £ 2 ^ and £ 2 < £ cr ~ 1.38 there 
are damping fluctuations near the minimum pq| . Let us note that in our case £ 2 < £ 2 r . 
To analyze the late time behavior one can linearize equation ( |3.1| ) as $ = $o — 5$ keeping 
only liner in 5& terms. A substitution yields the following equation for 5$ 

(-£ 2 d 2 + 1) e^d® = 35$. (3.2) 

At this point we want to mention two facts about this equation. First, if one omits the 
non-local exponential operator then one arrives to an equation for the harmonic oscillator 

(£ 2 «9 t 2 + 2)6^ = 

which does not give the desired late time behavior. Second, if we expand the non-local 
operator and keep only the constant and second derivative terms the situation depends 
on the value of £ 2 and for £ 2 > \ (which is our case) remains qualitatively the same. 
Using that e s * e~ mt = e m e~ mt one can try to find a solution to equation (3.2) of the form 
5$>(t) = Ce~ mt where for m we have the following characteristic equation 

(-£ 2 m 2 + 1) ei m2 = 3. (3.3) 

This is exactly equation (3.27) in |[^] with a redefinition m our = imp^j. 

Note, that this solution is valid only for late times. One expects thereby that the 
non-local operator plays a crucial role in the asymptotic regime. The latter characteristic 
equation for m is solved in Appendix [A[ Here only the results are used. A general solution 
for 5$ is an infinite sum 

= J2(A k e~ mkt + B k e mkt ) 

k 

because infinitely many roots m\ exist. This reflects the presence of an infinite number of 
derivatives. The latter means an infinite number of initial conditions which however are not 



5 See also |2^, |2^] for a cubic action. 

6 The potential is equivalent in terms of both <j) and <E> since (f> = $ for zero momentum. 
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arbitrary (see J2g] for a discussion on this point). Moreover in the full an action asymptotic 
behavior is subject to initial conditions in the origin. In other words considering the full 
non-linear equation and imposing initial conditions in the beginning of the evolution one 
arrives to a specific asymptotic configuration. We assume in our case a solution converging 
to an asymptotic value exists. For the definition of given in Appendix a vanishing 
condition imposes = for all k and a reality condition for 5& is formulated as = 
A*_ l _ k . Therefore, the most general real vanishing solution to equation ( |3.2| ) is 



6$(t) = J2( A k e ~ m " t+A * e ~ 



■m,t\ 



(3.4) 



fc>0 



Hence one yields a real function of time which oscillates with an exponentially decreasing 
amplitude. The function &{t) in turn goes to 1 oscillating with the decreasing amplitude 
near the asymptotic value. Recall that the solution constructed is valid only for late times. 
Let us also note that in a UV region, k > 1, one has to take into account string modes 
next to leading ones and this is beyond of our approximation. 

The main (i.e. the most slowly vanishing) contribution in (3.4) is given by k = and 
can be represented as 



S(p(t) = Ce sin(vt + <p) where 



1.1365, v ps 1.7051 



(3.5) 



1 rj2 

where we passed to 5 eft = es°« 5<F All other k will give a faster convergence and will play 
a role of small corrections. We find that a rather simple choice of constants C ~ 1 and 
ip ~ 2 in ( |3.5D matches the rolling solution presented in Fig. |l]a. Figures Fig. |b and Fig. |c 
contain corresponding plots. 




Figure 1: a. Solution to equation (3.1): $ (red line), 
b. An approximation to this solution 1 — S(f> given by i 
a different scale. 



(green line), <fi (blue line) for £ 2 = 0.9556; 
a) (blue line); c. The same plots as in b in 



Note that for £ 2 = (this case corresponds to a p-adic string 
equation (13.21) is simplified drastically and we have 



M with p = 3) 



4(log 3 + 2-Trki) 



where again different branches may be considered. The principal branch is k = and it 



corresponds to the rolling solution [25, 26, 27]. Other k do not match the rolling solution. 
They will give oscillations which were not present in a numeric analysis. Absence of such 
an oscillating behavior is related to initial conditions in the origin. One can also say that 
TTifc for k = 1,2,... are in a UV region beyond of our approximation. 
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The main conclusion of this Section is that the non-local operator e s * for positive 
s acts like a friction. For ^ 2 = the late time behavior of the tachyon is just a smooth 
rolling with a monotonically vanishing velocity to the minimum of the potential. This 
is not surprising mathematically but looks rather strange physically. Term with £ 2 in the 
l.h.s. of equation Q3,2| ) accelerates the tachyon but for actual £ 2 a friction becomes stronger 
and the tachyon eventually stops in the minimum. 

Another point of view is that for the considered £ 2 time intervals when the kinetic 
energy on space homogeneous configurations is negatively defined dominate providing a 
phantom behavior for the tachyon field. This justifies a phantom approximation used in 

B 



4. A Late time Rolling Tachyon in the FRW Universe 



We consider action ( p.2[ ) for the tachyon coupled to the gravity. We are going to employ 
an analog of the asymptotic expansion used above in Section [||. To this end we have to 
expand cf> = (f>o — 5<f> (and accordingly $ = (f>Q — 5<&, where <5$ = es v 5(f>) in Friedmann 
equations (|2.3j). In our notations <po = ±1 and the resulting equations read 



3H =7 { - 



e 

2 



6 2 + _5$2 + 6£i + §£2 + Aq 



8& 



(4.1a) 
(4.1b) 



where 



S£ 1 



dsd t e^ sV S^d t e-^ sV 5^ 



and Ao = — \ + A. The equation of motion for the tachyon field 
has the form 

(£ 2 £> + l)e--± v 5<$> = 3<5$. 



(2.4) after expansion 



(4.2) 



At this point we assume that in this approximation only the constant term Hq in an 
expansion of H survives in the latter equation. A validity of this assumption will become 
clear after a solution will be constructed. A value of Hq can be determined from equation 
( 4.1a ). However, even in this case an analytic solution to equation ( |4.2j ) in a closed form 
is not achievable so far. Instead of this we are looking for a solution constructed from 
eigenfunctions of the operator T>q = — 9 2 — SHodt- They have the following form 



5$ = Ae~ m+t + Be - *-* where m± = -H ± 



9H2 „ 
— - + m 2 . 



To become a solution to equation ( [4.2j ) with a constant H = Hq a parameter m should be 
determined by means of transcendental equation ( p.3| ) already appeared in the flat case. A 
general solution to equation (|4.2j) is 



5$ = ^ (A k e~ m + kt + B k e- m - kt ) 
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An analysis of this solution goes in an analogy with the flat case considered in Section 
For fh±k we have 

... _ 3 \0 k \ _ , v 
m±k = r± k + iv k , r± k = -Hq ± v k = sign(/3 fe )-, 



ml 



+ a k )+2 



'ml 



where a notation m\ = a k + if3k is used. However, in contrary with the flat case a selection 
of vanishing branches is not so obvious. It depends on a particular value of Hq. Notice, that 
v k is always non-zero, i.e. a solution is oscillating for any k. Fortunately, the symmetry 
m±k = persists. Thus a reality condition for 5& is A k = A*_ 1 _ k and B k = B*_ 1 _ k , 

and the most general real vanishing solution becomes 



5$ 



fc>0 



lite 



-m +k t + A* e -™* +k t 



+£'0 



k>0 



(4.3) 



where the prime means that only k for which the behavior is vanishing should be taken 
into account. 

Let us spend few lines discussing a question about the main (i.e. most slowly vanishing) 
contribution. A parameter Hq complicates the story since Hq > 0. Thus k = is not 
necessary the main contribution. Moreover, for specific values of Hq and k both plus and 
minus in f k may give a positive number. In this case a corresponding B k term is allowed 
and this is new compared to the flat case. To illustrate the situation rh k are drawn in 
Fig. for different values of Hq. On these plots each point corresponds to a specific k. 
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Figure 2: m±k for — (left), = i (middle), and Ho = 4 (right). Circles correspond to the 
plus and crosses to the minus sign. 

For k > the imaginary part is positive and grows with k. For k < — 1 the imaginary 
part is negative and grows in absolute value with |fc|. Circles correspond to the plus and 
crosses to the minus sign in ffi± k . Points with smallest positive real parts form the main 
contribution and one sees that a selection of such points is not so transparent. However, 
for the sequel we will need only to state that one can make such a selection. Thus for the 



main contribution in (4.3) we have 



Ce~ mt + C*e 



-m*t 



(4.4) 



- 8 - 



where we have passed to 6(ft = e s <5<J> and an index ±fc in m±k is omitted since a specific 
choice of k and a branch of the square root is made. Also for some values of Hq it is 
possible that two close points, say k and k + 1 for k > will have equal and smallest 
positive real parts. In this case they will only differ by a frequency of oscillations which 
will be suppressed with an equal exponential factor. However, for simplicity we will not 
consider this situation. Substituting solution (4^) into Friedmann equation ( 4.1b| ) one 
yields after an integration using relation (pi 



H = H + ^ (^V 2 ^ (1 + 4f - em 2 ) + ^fe~ 2 ^ (l + Af - fm* 2 ) ) (4.5) 



where an integration constant is chosen according to equation ( 4.1a| ) to be Hq — —■ < — 



do V 3 



and a constant C is the same as in (4^). Further one can check that equation ( 4.1a| ) 
is satisfied up to <50 4 terms which are beyond of our approximation. For the succeeding 
analysis it is convenient to represent (ft and H as follows 



1 - Ce~ H sm(ut + <p) and H = — J ^ - ^C H e~ 2H S m(2Pt + <p H ) (4.6) 

9o V 3 gi 



where C and <p are arbitrary (integration) constants, 



C 



■2 



C H = — Vmm*J(l + 4£ 2 ) 2 + £ 4 m 2 m* 2 - £ 2 (1 + 4£ 2 )(m 2 + m* 2 ), 
32 * 



_ 7r / m — m \ ( 

and ipu = 2ip -\ h arctan i— — — arctan % 

2 \ m + m* J \ 



ir f.m — m*\ (. ^ 2 (m 2 — m* 2 



2 + 8£ 2 - £ 2 (m 2 + m 



* 2 x 



In a special case £ 2 = which corresponds to a p-adic string in the FRW Universe the 
story is much simpler because eigenvalues fh±k can be expressed in elementary functions 
since m\ = 4(log3 + 2irki). Moreover, fh±o are pure real and m+o corresponds to a 
smooth approach of the tachyon field (ft to the asymptotic value (fto = 1 without oscillations. 
Expressions for 5 (ft and H in this case are as follows 

= Ce~ M and H = -\f^ - ^C H e- 2m (4.7) 
9o V 3 gi 



where m = |i?o + y^f 11 + 4 log 3, Ch = ^j<p and Hq is the same as in (fl~5|). Again one 



can check that equation ( 4.1a ) is satisfied up to <5</> 4 terms. 

Plots representing the tachyon field evolution both for £ 2 ~ 0.9556 and £ 2 = are 
shown in Fig. [||. 

5. Cosmological consequences and further directions 

The obtained asymptotic solutions for the tachyon and Hubble parameter ( |4.6| ) lead to a 
number of interesting cosmological properties. 
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First, we mention that H tells 
us about the time point to after which 
the solution is reasonable. Indeed, 
for H to be strictly positive we have 

to require */4p > ^C H e- 2H . Solv- 
ing this inequality w.r.t. the time we 

havei>t = -^log(W 
This is the characteristic time in ques- 
tion. 

Then one sees that during the 
late time evolution of the tachyon 

field the Hubble parameter goes to a constant. H obviously vanishes and moreover, both 
H and H do oscillate. The state parameter w also has an oscillating time behavior and 
goes asymptotically to —1. An explicit expression for it is of the form 




Figure 3: Tachyon field for £ 2 « 0.9556 (left) and £ 2 = 
(right) for C = 1, (p = and A = 3 • 1CT 4 . 



-2r sin(2i4 + <£>#) + 2v cos(2^t + </?#) 



f o C H e- 2ft sm(2vt + (p H ] 



(5.1) 



It is very interesting that w crosses the phantom divide w = — 1 during the evolution. Such 
a crossing is forbidden in single field cosmological models with a local action. In our model 
we see that a non-locality breaks this restriction. The deceleration parameter q behaves 
very similar to w because its expression through H is very close to w 



q = -l 



-2r sin(2z^ + (pn) + 2v cos(2i4 + ipu) 
^ - ^C H e- 2ft sm(2ut + <p H 

•J yo x 



(5.2) 



Hence the Universe exhibits an acceleration but because of oscillations quintessence and 
phantom phases change one each other with the time. 

The scale factor a is related to H as H = a /a and can be readily found to be 



a e 



J Hdt 



clq exp 




n - 2 ft rsi ' a { 2ut + L PH) + vcos{2vt + ipH)\ , K „, 
C-^e - 2 , - 2 (5.3) 



where ao is an integration constant. For late times it has a simple exponential approxima- 
tion 



a = a^eso v 3 _ 



To plot H and w we specify Hq (or Ao) to be so small that the f and v in fl4.6| ) are 
f+o and P+o respectively. This corresponds to the middle plot in Fig. [|. The explained 
above behavior of H, w, q, and a is visualized in Fig. ^. We point out that in spite of 
the presence of only one scalar field we observe similarities in the behavior of cosmological 
quantities with the two-field model analyzed in [31]. These are crossing of the phantom 
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Figure 4: From left to right H, w, q, and a for n = g = C = do = 1, (p = and Ao = 3 • 10 21 . 

divide by w and the qualitative form of plots for H, uu, and q. An order of magnitude of 
Ao comes from an analysis performed in Ej| . 

For £ 2 = one readily gets using (|4.7| ) for the state and deceleration parameters 

2C H me- 2fht 2C H me- 2rht 



w 



■1 



-2mt 



9o 



2 • 



1 



4a - iiC^e- 2 ™* 

■J 5o 



2 • 



(5.4) 



We see that w is always less then — 1 in this case and a crossing of the phantom divide does 
not occur. Thus the tachyon behaves purely like a phantom for £ 2 = 0. The scale factor a 
in this case is found to be 



a = ao eJ Wd * = a exp — W-^t + — -4e 2mt » a efoV 3 



(5.5) 



where ao is an integration constant. Plots for H, w, q, and a in the case £ 2 = are 
presented in Fig. ||| The behavior of the model in this case is very close to a single 





Figure 5: From left to right H, w, q, and a for n — g = C = ao = 1 and Ao = 3 • 10 21 in the 
case £ 2 = 0. 

phantom field model solved in [Hi]]. 

Remarkably, that both for £ 2 ~ 0.9556 and £ 2 = the Big Rip singularity problem 
is avoided because w exhibits a non-trivial time dependence and consequently is not a 
constant less then —1. 

Constructed solution reveals many properties expected from the cosmological point of 
view giving rise to an interest to a further investigation of this model as well as models 
coming from the SFT in general. Apart from many other possible directions of developing 
this type of models we would like to mention ones which are seen as the most important. 
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It would be very interesting and at the same time very difficult to construct a solution to 
full non- linear equations (|2.3|). It seems that if it can be done it will be a numeric solution. 
However, even a numeric approach faces a difficulty of dealing with an exponential of a 
Beltrami-Laplace operator in a curved background. The main technical problem is that a 
Beltrami-Laplace operator itself contains an unknown function H to be determined while 
solving the equations. 

A complete analysis of a stability of the obtained solution is also of great importance. 
A related issue which is an inclusion of other cosmic fluids and especially the Cold Dark 
Matter (CDM) forming about 1/3 of the present day Universe and an investigation of a 
dynamics of such a coupled model and its stability would be very important (see [42] as 
an example of a coupling of one phantom field to the CDM). 

Also it would be interesting to find an exact solution to full equations with infinitely 
many derivatives in curved background probably by adding extra terms to the action. 
This is in accord with an analysis performed in JO], 43] where a small in terms of coupling 
constants correction to the potential makes the problem analytically solvable. An existence 
of an analytic solution provides a possibility of a qualitative analysis without a strong 
support of numeric methods. 

Another important question is a consideration of closed string scalar fields, the closed 
string tachyon and dilaton, as well as their coupling to the open string tachyon. This 
problem is related to finding of lump solutions and a development in the flat background 
was started in [43|. Similar lump solutions have been constructed in 44]. An extension to 
the FRW Universe using a linearization of equations of motion developed in the present 
paper would be an interesting analysis of the role of closed string excitations 7 . 

In spirit of a selected role of vector fields in a construction of local phantom models 
without an UV pathology |24|] it would be very interesting to incorporate a string vector 
field that has in the SFT approach a non-local interaction in a study of a nonlocal rolling 
tachyon dynamics. 
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A. A Solution to the characteristic equation 

The equation in question is the following 

7 See [M for a discussion on the closed string tachyon and dilaton condensation. 
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It is a transcendental one and a solution to it can be represented analytically as 
m l = ^ + 4W k (-4^2^) = 4A + m * (-3Ae- A ) 

where W k is the k-s branch of Lambert W function satisfying a relation W(x)e w ^ = x. 
It has infinitely many branches distinguished as W k . A branch analytic at is referred 
to as a principal one. There is exactly one such a branch and k = is assigned to it. 
The branch cut dividing Wo and W±\ is the part of the real axis (— oo, — 1/e) and — 1/e 
is the branch point. The branch cut dividing all other neighbor branches is the negative 
real semi-axis and is the branch point. For the principal branch Wo(x) the image of the 
(— oo, —1/e) interval is — /3cot (3+i{3, where (3 G (0, ir]. For the branch W-i(x) the image of 
the (— oo, —1/e) interval is — /3cot (3 + i(3, where (3 E (— it, 0]. For all other branches W k (x) 
the image of the (— oo, 0) interval is — /?cot (3 + i(3, where (3 6 (2kn, (2k + l)ir] if k > and 
(3 G ((2k + l)-7r, (2k + 2)-7r] if k < —1. For our purposes we are interested in the argument 
x = — 3Ae _A ~ —0.6042 < —1/e, so the above properties are relevant. To understand the 
dependence — (3 cot (3 + i(3 we are going to consider the equation 

ye y = x 

for x < —1/e. Assuming a complex solution y = a + i{3 we have 

x = ( a + i(3) e a+il3 = ^a 2 +p2 e a e< aicttm ° + P). 

Since x is a negative real number we have to require a relation arctan ^ + (3 = (2k + l)ir 
which yields a = — /3cot (3 where a branch of the cotangent is the origin of branches of W 
function. Using this expression for a we can write x = — ^j^ e ~' 3cot/3 - Here the dependence 
on (3 is manifestly odd. From the point of view of W function the change (3 — ► —(3 is the 
change of branches k — > —1 — A;. This means the relations 

ReW fc (x) = ReW_i_ fc (x), lmW k (x) = -ImW_i_ fe (x) 

hold for x < —1/e. Moreover absolute values of real and imaginary parts of Wk(x) grow 
when k grows for k > or k decreases for k < —1. The latter is not evident and can be 
derived from the series expansion of W function. 

For our value of £ 2 all roots are complex with non-zero real and imaginary parts. 
Denoting m\ = au + i(3k we have 



m k = r k + iv k , r k = —— , v k = sign(/3 fe )-u, v = J -2a k + 2^ a\ + /?| 

where for all square roots the arithmetic branch is chosen. In this case a k > 0. Another 
root is m k = —(r k + iv k ). Using the symmetry a k = a-\- k , (3 k = —(3-\- k we have for our 
definition of m k that m k = m*_ 1 _ k . Also for our value of £ 2 all a k < 0. 
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